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• Abstract 

' There are many evolution partial differential equations which can be 
cast into Hamiltonian form. Conservation laws of these equations are re- 
lated to one-parameter Hamiltonian symmetries admitted by the PDEs [2- 

^ , The same result holds for semidiscrete Hamiltonian equations 2 . In this 

^\ ' paper we consider semidiscrete canonical Hamiltonian equations. Using 

. symmetries, we find conservation laws for the semidiscretized nonlinear 

' wave equation and Schrodinger equation. 

o , 

^ ■ 1 Introduction 
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Many partial differential equations (PDEs) of nondissipative continuum me- 
chanics can be presented in Hamiltonian form (see and references therein). 
It is well known that conservation laws of Hamiltonian PDEs are related to 
one-parameter Hamiltonian symmetries The analog of this result holds in 
the semidiscrete case with no change in the statement although the framework 
must be modified 

^ ] In the present paper we examine canonical Hamiltonian equations 

which form a special type of the Hamiltonian equations. For instance, the 
nonlinear wave equation 

vtt = Av- V\v), (1.2) 
and the nonlinear Schrodinger equation 

#t + AV' + F'(|^P)V = 0. (1.3) 

can be rewritten in the form We will introduce semidiscrete analogs of 

the canonical Hamiltonian PDEs and show how one can use the Hamiltonian 
form of Noether's theorem to find conservation laws of these equations. 

To find conservation laws of semidiscrete equations with the help of Noether's 
theorem, we need to know symmetries in the evolutionary form. Different meth- 
ods which can be used to find symmetries of discrete and semidiscrete equations 
are discussed, for example, in The most successful application was shown 
for linear difference equations where for the considered equations 
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symmetry algebras isomorphic to those of the underlying continuous equations 
were found. However, it is not known how to find all symmetries admitted by 
nonlinear semidiscrete equations. 

One of the possibilities, which we exploit in this paper, is to make use of 
the admitted Lie point symmetries. Many Lie point symmetries of semidiscrete 
equations can be easily found as Lie point symmetries of the underlying con- 
tinuous equations preserved under the space discretization. In both continuous 
and semidiscrete cases these symmetries are given by the same vector fields. 
Using factorization, one can obtain corresponding evolutionary operators. If 
the latter are Hamiltonian symmetries, they let us find conservation laws. Ob- 
viously, for our purpose we will not be interested in all symmetries but only in 
Hamiltonian ones. 

The layout of the paper is as follows: In Section 2 we briefly introduce Hamil- 
tonian equations and specify canonical Hamiltonian equations. Symmetries and 
the Hamiltonian form of Noether's theorem are discussed in Section 3. In Sec- 
tions 4 and 5 we examine semidicretizations of equations (|1.2|) and (|1.3|) and 
find their conservation laws. In final Section 6 we make concluding remarks. In 
particular, we mention the connection between Euler-Lagrange equations and 
canonical Hamiltonian equations. 



2 Hamiltonian equations 

For simplicity we will consider the case of one space coordinate x. We assume 
that the solutions are sufficiently smooth, all variational derivatives tend to 
zero as the solution tends to zero and the solution and a number of its space 
derivatives tend to zero as |x| — > oo. We suppose that the solution decays fast 
enough so that all integrals and sums make sense. 



2.1 Hamiltonian partial differential equations 

Many systems of evolution equations 

= i^(x,u(-)), 

where u stands for dependent variables u = (u^ , ...,u^)'^ and u^"^^ = 
(u, ui, U2, Urn) represents u and a finite set of derivatives of u with respect 
to space coordinate x, can be cast into the Hamiltonian form 

(2.1) 

with the help of the Hamiltonian functional 'H[u], variational operator 5 ■ /6u 
and the linear operator V p. We denote as T the space of functional 



The operator V must be Hamiltonian, i.e., it forms the Poisson bracket 
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satisfying the conditions of skew- symmetry 

{V,C] = -{C,V] (2.3) 

and the Jacobi identity 

{{V, £}, 7^} + {{7^, v},c} + {{/:, n}, v] = o (2.4) 

for all functionals V, C,Tl £ T . 

The variational derivatives of a functional can be found by the action of the 
Euler operators on the integrand 

= i=l,...,]V, 

where Dx is the total space derivative operator. 



2.2 Canonical Hamiltonian equations 

Canonical Hamiltonian equations form a subset of the equations char- 
acterized by an even-dimensional space of dependent variables = 2n, u = 
(u^, f", if^, ...,t(;")"^, and the canonical Hamiltonian operator 

■^=(-'»o;;)- 

where /„ is the n x n identity matrix and 0„ is the n x n zero matrix. Thus these 
Hamiltonian equations have the form It is easy to see that the Poisson 

bracket generated by operator J 



satisfies skew-symmetry ()2.3() and Jacobi identity (|2.4() . 



2.3 Semidiscrete Hamiltonian equations 

Given a Hamiltonian PDE, we attempt to discretize both the Poisson bracket 
and the Hamiltonian functional so that we preserve Hamiltonian structure. 

To consider semidiscrete equations we introduce a two-dimensional mesh 
which is uniform in space and continuous in time. Let us denote the mesh 
points as {xi{t)}, i £ Z, t > and define mesh Q by two conditions: 

U : Xi+i{t) - Xi{t) = Xi{t) - Xi-i{t), Xi{t + T) = Xi{t), i£Z, t,T>0. 

(2.7) 

The first equation requires the space mesh to be uniform for any fixed time. 
The second equation requires that only vertical mesh lines in the time-space 
plane are considered. 
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Now we can introduce discrete space derivatives u\ = j)u^ ^ u\ = Du\, 

h +h h -hh 

"2fc+i = DuL, linfc+o = DuL i = I,..., I, where D and D are the right 

h ^ +hh h ^ -hh ^ +h -h 

and left discrete differentiation operators 

D = ^j—; D = , 2.8 

+h h -h h 

defined with the help of the right shift Sj^ and left shift operators 

S+f{x) = f{x + h), S.f{x) = f{x-h). (2.9) 
We will consider the space of discrete derivatives u^"^) = (u, u i, u 2, u 

h h h h 

and functionals of the form 

p = yp(t,x,/i,u(-))/i, (2.10) 

where the summation is taken over all space points of mesh for some fixed 

time. We denote the space of such functionals as !F. 

h 

We assume that the Hamiltonian operator T) can be approximated by an 
operator V such that the discrete bracket 

h 



n 



defines a Poisson bracket for functionals from J^, i.e. the bracket {•, ■}h is skew- 

h 

symmetric and satisfies the Jacobi identity. Then we can choose some approx- 
imation fi of the Hamiltonian 7i to obtain the set of semidiscrete Hamiltonian 

h 

equations 

which approximate equation ()2.1|) on mesh ()2.7() . 

For the space discretization of the canonical Hamiltonian PDEs we can keep 
the canonical operator J since it generates a discrete Poisson bracket, namely 



{V,C} 



h 




5V6C dvSC" 

h h h h 



h h ^ \ 6v^ 6w^ 5w^ 5v^ 



and take a discretization fi of the Hamiltonian functional Ti.. This procedure 

h 

provides us semidiscrete canonical Hamiltonian equations 

sn 5n 

V, = -^, = --A, jeZ, (2.13) 

where we use vector notations v = {v^, ...,v^)'^ , w = {w^ , ...,w^Y' . 
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3 Symmetries and conservation laws 
3.1 Invariance of semidiscrete equations 

Let Z be the space of sequences of variables (t, x, /i, u, u i, U2, ...) and A be the 

h h h h 

space of analytic functions of a finite number of variables z from Z- 

h 

Invariance of the semidiscrete equations 

u = F{z), Fie A (3.1) 

h 

defined in the points of some two-dimensional mesh was considered in 
Symmetries of equations (|3.1j) are transformations generated by vector fields of 
the form 

X = eiz)-^^ + e{z)^ + ri\z)-^^ + • • • , ^-,rj^ G A, (3.2) 

which leave the equations and the mesh invariant. The infinitesimal criterion 
of invariance can be presented by the three conditions 

X{u-F{z)) =0; (3.3) 



DD{n=0; A(n=0; (3.4) 

which are to be satisfied on the solutions of Condition ()3.3() requires the 

invariance of Eqs. (|3.1j) while conditions (|3.4jl that of mesh il.. The operator X 
must be prolonged on all variables appearing in the equation p.lf) 

X = - + S,^ ^ + ff —— + (jf—— + C * ^ + C * ^ -\- • ■ ■ -\- D (3 5) 

dt dx du^ dv} du\ u"^ du\ +h dh 

On the uniform in space grid the coefficients of the prolonged operator are 
defined by the prolongation formulas: 

cP' = Dtirj') - iiWtie) - niA(e), Ci = DW) - S+{u')D{e) - n\D{e). 

h +h +h h +h 

ci = DD{r]') - 2 uiDie) - }iS+{u')D{e) + j^S.{u')D{e), ■ ■ ■ 

h —h+h ri +h +h —h 

(3.6) 

Note that u\ = Dx{u^) is the "continuous" derivative. It is supposed to be 
in some discrete representation, for example, Dq{u), which will be introduced 

h 

below. 

Operators of the form (|3.2|) are called Lie-Backlund (or generalized) symme- 
tries. It is a difficult task to find Lie-Backlund operators admitted by discrete 
equations. However, Lie point symmetries 

Ci C\ C) 

X = e{t,x,u)-+e{t,x,u)—+l^\t,x,u)^ + ■■■ , (3.7) 
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where coefficients , and r] depend only on dependent and independent 
variables, are easier to detect since such symmetries are given by the same vector 
fields in the continuous and discrete cases (7| • Practically, one can check whether 
Lie point symmetries admitted by the underlying PDEs are admitted by the 
semidiscrete equations or not. Although this procedure does not guarantee that 
we find all symmetries, it let us avoid solving discrete determining equations. 



3.2 Factorization of operators 

Following |7j,[S], let us consider a special operation of left multiplication of a 
Lie-Backlund operator by an analytic function ^(z) G A- 

h 

e * ^ = k'-^^ + ie^ + ev/i + • • • + (Ir )|-. (3.8) 

ot ox Ou^ +h oh 

The first coordinates in operator are multiplied by ^{z) while the remaining 
coordinates are computed according to the prolongation formulas (|3.6() . 
The operator 

i\z)*Dt + e{z)*D^, (3.9) 

h 

where Dt is the total time derivative operator and D x is a discrete presentation 

h 

of the total space derivative operator Dx^ plays the role of an ideal of the 
Lie algebra of operators [2]. There are several possibilities to choose 

operator D x- One can take a representation based on the right or left discrete 

h 

derivative |ZI,[H] 

rln' I u rill I u ^— • n i u 



(3.10) 



D- = — + D{u)— + ..., D = ^ D", 

ox -h ou -h ^ n -h 

n=l 



or the discrete representation based on the central difference derivative |U 

A:=0 

(3.11) 



with coefficients 

, ,.13 2k-l 1 (-1)'= (2A; - 1)!! . ,„ . . 

a2k+i = (-1)'' : : = \ ^ TTT^, 2A:-1 !! = 1-3-5 • • • (2k-l). 

Let us mention that operator Dq can be presented in terms of powers of 

h 

the shift operator 

oo 

^0= CkSk, Sk = {S+)'' (3.12) 

h ^ — ' 
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with convergent coefficients Ck- To find this representation we rewrite this 
operator as 

^ oo 

?«=2^E^(«.-S-)--. (3.13) 

fc=0 

Using 



lfc+l+p. 



(2A; + 1)! 



xfc-p (2A: + 1)! 
p=o + '^^'^^^ ' - p)l{k + 1 + p)l 



5'-2p-l 



we obtain 



1 °° 

I* = ^ X] Cp(5'2p+i - S-2p-i), (3.14) 



p=0 

where the coefficients are given by the series 

oo 



r - (-UP\^(-^)k^2k+l^_{2k + l)]_ 



^A: + l+p22fc(A;-p)!(A;+p)! 



(3.15) 



Lemma 3.1 T/ie series icl.l5\) defining coefficients Cp forp £ N are convergent. 
Proof: Using (A; + p)\{k — p)l > (kl)'^ we get 

2 



Fpl < 



£:^fc + i+p V(2'=(^)0V 



fc=p 

We use the bounds on the factorial provided by the Stirhng's expansion 

n" exp(— n)V^27m < n! < n"\/27rnexp ( — n H ) for nGN 

\ 12n / 

to obtain the inequahty 



< ——: — r exp ( — ^ ) for G N 



1 (jmLr<—^ -..J' 



k + l+p\{2''{k)iy J TTk{k + l+p) \l2k^ 

that insures us that the series defining Cp converges, moreover Cp ~ l/|p| as 
p ^ oo. 

□ 

One can easily check that the operators D and D can not be presented in 

+h -h 

the form (IIH2I) with convergent coefficients. 
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Lemma 3.2 The operator D q is skew-adjoint , i.e. 

h 



ViDoUih = - y^ UiDoVih 



h ^ h 

j=— oo j=— oo 

for {ui}, {vi}, i €z such that Ui,Vi ^ as i ^ oo . 

Proof: The equality can be check with the help of (|H.14j) . □ 
In what follows we will use D^. With the help of the ideal * Dt + * 



we can find the evolutionary operator corresponding to the operator (|3.2|) 

X = X-e**A + r*^° = ¥T?- + --- , fj' =v'-^^u'-C''Doiu') (3.16) 

au* h 

On the solutions of equations H2.12() we can exclude time derivatives u'^, i = 
1, ...,N from the coefficients of the evolutionary operator. 

It is important to note that although the factorization let us find evolu- 
tionary operators corresponding to given Lie point operators, generally, the 
obtained evolutionary operators does not have to be admitted by the semidis- 
crete equations which admit the original nonevolutionary operators. In the 
continuous case ^* * = • Dt and .^^ * = S,^ • D^, i.e., the operation * 
is equivalent to a left multiplication of the prolonged operator 0. The oper- 
ators Dt and D^ are admitted by all differential equations and, consequently, 
the operators ^* * Dt and * Dx are also admitted by all equations. It fol- 
lows that an evolutionary operator is admitted if it corresponds to an admitted 
nonevolutionary operator p. 

This result does not hold in the semidiscrete case. In the general case 
operator * Dt = ^* • Dt is admitted, but operator * D^ is not. Thus, having 
obtained an evolutionary operator, one has to check that this operator or used 
in the factorization operator * D^ is admitted by the considered semidiscrete 
equations. 

We can provide only a very restricted class of operators * D^ which are 
admitted by an arbitrary discrete equation. 

Lemma 3.3 For such that D{i^) = we have 

+h 

*D° = ■ D° 

and, consequently, this operator is admitted by all difference equations. 
Proof: The result follows from the prolongation formulas. □ 
Remark The condition D (C^) = is very restrictive, but its multi-dimensional 

+h 

analog 

D = 0, 

where Xi is an independent space variable and D is the right discrete derivative 

+hi 

with respect to Xi, leaves more freedom. For example, it allows rotations in 
XiXj , i ^ j planes 

X - A_ A 

■^i.j — Xi Xj 

OXi OXi 
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3.3 Conservation laws 

For a system of Hamiltonian equations (|2.13jl considered on the grid ()2.7j) we 
have the fohowing types of the conservation laws: 

1. Conservation of symplecticity. 

Due to their canonical form semidiscrete equations possess the conservation 
of symplecticity 

^ oo oo n 

—oj = 0, w = ^ dvj A dwj h= ^ '^dvj A dw] h, (3.17) 

j=—oo j=—oo i=l 

where c?Vj = {dvj , ...,dVj')'^ and dwj = {dWj, ...,dwj')'^ are solutions of the 
variational equations 





fc=i 

We suppose that dvj, dwj ^ as j — > cxd. Differentiating the 2-form w 

^ oo n oo n 

—uj = dv^, A dw] h+ "S^ dv\ A dw] h 

j=— oo i=l j=— oo 1=1 

= y y du;tiAd< h+ y y \ 1 dv'^.^iAdw', h 

j,(=— ooi,fc=l \ J / j,(=— oo i,fc=l 

j,t=— oo j,fc=l J+' \ J I j,l=—ooi,k=l J+' \ J 





d 



h \ h \ O h 



dWj_^_i \ Sw'j j dwj \ Swj^i j ' 9vj_^_i \ 6vj j dvj \ ^v'j^i j ' 
d 



6n 

h_ 




we obtain the conservation of the symplectic form. 

This is a generalization of the symplectic structure for the Hamiltonian 
ODEs 9 to the infinite set of semidiscrete equations. Let us note that recently 
proposed multi-symplectic formulation of PDEs jn]-^ni allows to consider 
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local conservation of symplecticity, which we do not have in the present frame- 
work. 

Underlying equations possess the conservation of symplecticity 



dt 



LU = 0, iu = J dv Adw dx = j '^dv'- A dw'' dx (3.18) 

i=l 



that is the continuous limit of H3.17|l . In this case dv and dw are solutions of 
the variational equations for 



k=l I L ^ ' k=l I ' 

fe=l I I ^ ' ; t 



fe=l / t ^ ^ fc=x ; J ^ ' 

where the operator d denotes the vertical differential of the vertical form it 
acts on. In our case we have vertical 0-forms, i.e. functions. We assume that 
dv, dw and their space derivatives which appear in the variational equations 
decay at a; — > (X). Let us show that the conservation of symplecticity using the 
variational complex (see, for example, ^): 

—uj = / f ^ d-Uj A dw'' + ^ dv' A dwl \ dx 

\i=l i=l J 



-^dw^ I dx 



Using the integration by part, we obtain that 

/ E "B-^-' = / E E(-".)'iJ*' = / E E If *i 

i=l i=l k ^ i=l k ^ 

because du^, i = 1, ...,n tend to zero as x ^ oo. Since the same result is valid 
for variables w we get 



dH dx = 6n 



so that 

= s'n = o 

dt 

as follows from the exactness of the variational complex. 
2. Conservation of distinguished functionals 
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Definition For a given Hamiltonian operator T> a distinguished functional is a 

h 

functional C{x, h, u^")) such that 

h h 



(3.19) 



It follows that a functional is distinguished if and only if its Poisson bracket 
with every other functional is trivial: 

{C,n}h = for any H ■ (3.20) 

h h h h 

Distinguished functionals are conserved by semidiscrete equations originat- 
ing from any Hamiltonian functional. For the canonical bracket there are no 
nontrivial distinguished functionals so that we do not get conservation laws of 
this type. 

3. Hamiltonian form of Noether's theorem. 

Definition The Hamiltonian vector field associated with a functional "P is the 

h 

unique smooth vector field Xp satisfying 

Xp{j^) = {j^,V}h (3.21) 

h h h 

In the coordinate form it can be presented as the operator 

h 



Some symmetries ()3.2() are given as Hamiltonian vector fields or are equiv- 
alent to Hamiltonian vector fields under factorization (|3.16j) . Such vector fields 
let us use the following theorem 

Theorem 3.4 For a Hamiltonian system of semidiscrete evolution equations 

112.1^) a Hamiltonian vector field Xp determines a generalized symmetry of the 

system if and only if there is an equivalent functional V = V — C, differing from 

h h h 

V only by a time-dependent distinguished functional C{t,x,h,u^^^), such that 
h h 

V determines a conservation law. 

h 

For the operator J generating the canonical bracket a Hamiltonian vector 
field has the form 

where V is the generating functional. The canonical bracket has only trivial 

h 

time-dependent functionals, i.e. functions f{t). Thus, for a Hamiltonian sym- 
metry (I3.23|) there corresponds a conservation law V = 'P—f{t), where the 

h h 
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function f(t) needs to be found with the help of the considered equations. In a 
particular case when the Hamiltonian functional and the considered Hamilto- 
nian symmetry are time independent we obtain a linear function f(t) = at + b, 
a,b = const. 



4 The nonlinear wave equation 

In this section we consider the nonlinear wave equation 

vtt = - V'{v), (4.1) 

where V{v) is some smooth function. For simplicity, we consider the case of 
scalar v. With the help of a new variable w = vt, the equation (|4.1() can be 
rewritten as the system 

vt = w; 

(4.2) 

Wt = Vxx - V'{v). 

This is a canonical Hamiltonian system generated by the Hamiltonian functional 

^ = / (t + 1 + ^^''O ^^'^^ 

Let us take the following approximation of the Hamiltonian functional 

nh = Y.H[v,w]h, H = ^ + ^ + V{v). (4.4) 

Q 

It provides us the system of the semidiscrete equations 



w = V 2 — V'iv). 



(4.5) 



h 



For arbitrary V{v) the admitted transformation group for (|4.5|) is two- 
dimensional. Its Lie algebra is spanned by the operators 

d d 

Lorentz transformation 

Xs = x— + t—, 
at ox 

admitted by underlying system (|4.2() with arbitrary V{v) is lost under 

discretization since it brakes the mesh invariance. Thus in the case of arbitrary 
V{v) Noether's theorem gives two conservation laws: 

1. The time translation Xi leads to the conservation of the Hamiltonian 
functional fi. Indeed, the factorization of operator Xi gives the evolutionary 

h 

vector field 

X^ = nj^+(v2-V'iv)^ ^ 



dv \h I dw ' 



12 



which is generated by the Hamiltonian functional. Thus we obtain conservation 
of Hamiltonian functional (|4.4() that stands for the conservation of energy. 

2. The space translation X2 corresponds to the evolutionary operator 

Checking condition 1)3. 23(1 for the coefficients of operator X2, we find the gen- 
erating functional 



V2 = ^wDQ{v)h = -y2'^^o{w)h. (4. 



This functional is a conservation law of the semidiscrete equations 1)4. 5|) . In the 
continuous limit it corresponds to the functional 

V2 = j wvxdx = — j vwxdx, (4.9) 

which is a conservation law of the Eqs. ()4.2|) . In physical applications this 
conservation law is referred as linear momentum. 

In order to obtain this conservation law we need to consider "nonlocal" 
functionals, i.e. functionals which are defined as double sums over the mesh 
points, since operator D q is a sum over infinitely many mesh points. 

h 

Let us consider the special cases of the potential V{v) which lead to addi- 
tional Hamiltonian symmetries 

2 

a) For the quadratic potential V{v) = there is an infinite series of con- 
servation laws for the wave system ()4.5() . In this case the semidiscrete system 
admits the infinite set of Hamiltonian operators 

d d 
Yk = Do{v2k)^ + Do{w2k)^, A; = 0,1,..., 

h h OV h h OW 

which provides us the infinite set of the conserved functionals 
Rk = yZ'wDo{v2k)h = - yZvDoiw2k)h 

In the continuous limit these conservation laws correspond to the functionals 

Rk = j VW(2k+l)dx, 

which are conserved quantities for the system 1)4.2(1 . 

Since for the quadratic potential we get the system of linear equations their 
solutions possess a superposition principle. It is reflected in the invariance with 
respect to the symmetry 

Z = a{t,x)^ + at{t,x)^, 
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where the function a{t, x) is an arbitrary solution of the equation 

, , a(t,x + h) — 2a(t,x) + a(t,x — h) ^ , , /.-^n 
atf (t, x) = ^ ^ |-f^ ^ - Ca{t, x). (4.10) 

The operator Z is Hamiltonian. It corresponds to the functional 

T=/ (a(t,x)w — at(t,x)v)h, 

which is a conservation law of the semidiscrete linear system if function a(t, x) 
satisfies the equation H4.10() . In the continuous limit this functional goes into 
the functional 

T= {a{t,x)w — at{t,x)v)dx, 



which is a conservation law of ()4.2() if the function a{t, x) satisfies the equation 

att{t, x) = axxit, x) - Ca{t, x). (4.11) 

b) Let us consider the case C = 0, i.e. the wave system (|4.5|) without a 
potential V'{v) = 0, in details. The symmetry Z is specified by a solution of 
the equation 

^ _ a{t, x + h) - 2a{t, x) + a{t, x - h) 

We present a number of the conservation laws corresponding the symmetry Z 
taking particular solutions a{t,x) of the equation 1)4.12(1 in Table 2. 

Let us remark that the continuous conservation laws in the case V'{v) = 
are given by the function a{t, x) satisfying the equation ()4.11|) with C = 0. The 
general solution can be written down as 

a{t, x) = ai{t — x) + a2{t + x), 

where ai and a2 are arbitrary functions. 



Example 4.1 As we mentioned before the factorized operator may fail to be 
admitted. Let us consider V'{v) = 0. In this case system (|4.5() admits the 
scaling symmetry 

d d d 
= -t— - X— + w—. 
ot ox ow 

The corresponding evolutionary vector field 

d ~ d 

X^: = {tw + xDo{v))— + {w + tv2 + xD o{w)) — , 

h OV h h OW 

found with the help of factorization formula (|H.lHj) . is Hamiltonian. It is gen- 
erated by the functional 



V* =tn + y^ xwDo{v)h. 
h ^ n ^ 
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However, neither the functional P * is a conservation law nor the canonical 

h 

operator X^, is a symmetry of the semidiscrete system (j4.5j) with V'{v) = 0. It 
happens because the operator 

ox h OV h OW 

which is used for the factorization 

X^ = X^-t*Dt-x*D^, 

is not admitted by the semidiscrete system: the second equation does not allow 
this operator since 

2v2 + xDo{w) ^ DD{xDo{v)) (4.13) 

h h -h+h h 

on the solutions of the semidiscrete equations. In the limit /i — > ()4.13|) turns 

into an equality and the continuous limit of P * , namely the functional 

h 



V* = t Ti + J xwVxdx, 
is a conservation law of the system (|4.2|) with V' = 0. 







5 The nonlinear Schrodinger equation 

Another equation which can be cast into the canonical Hamiltonian form is 
the nonlinear Schrodinger equation ()1.3() which arises in nonlinear optics. It 
describes the main features of the beam interaction with a nonlinear medium 
and is considered as the basic equation of the nonlinear optics |15j . The equation 
also has important applications in plasma physics jl6j . 
Let us consider the case of one-dimensional space 

ii^t + i^xx + F'{\^p\^)i; = 0. (5.1) 

For real and imaginary components v and w {ip = v + iw) the Schrodinger 
equation can be rewritten as the system 

Vt = -Wxx - F'{v^ + 'U?)W] 

(5.2) 

Wt = Vxx + F'{v'^ + w'^)v; 

which is a canonical Hamiltonian system. It is generated by the Hamiltonian 
functional 

n = U {\i^x\'-FM'))dx. (5.3) 
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The system of equations (|5.2() in the case of arbitrary F admits a four- 
dimensional transformation group presented by the operators |15j : 

d d d d d d d 

Xi = —, A2 = — , X3 = w- V—, Xi = 2t- XW— + XV—. 

at ox ov aw ox ov ow 

(5.4) 

2 

In the case F[z) = there is an additional scaling symmetry 

^ d d d d . ^. 

X5 = 2t— + x—-v—-w—. 5.5 

ot ox ov ow 

Physically, this particular case of F{z) corresponds to an isotropic medium with 
the cubic polarizability (first approximation of the nonlinear polarizability) . 
Let us discretize the Hamiltonian functional as 

then we obtain the system of semidiscrete equations 

i) = — W2 — F'{v'^ + w^)w] 



W = V2 + F'{v'^ + v?)v. 



(5.7) 



Now we are in a position to go through the symmetries (|5.4j) . (|5.5|) to check 
whether they are preserved under the space discretization. If they are admitted 
by the semidiscrete system (j5.7|) , they may provide conservation laws according 
to Theorem 13.41 

1. The time translation Xi is admitted by system (|5.7|) . The corresponding 
to Xi evolutionary symmetry is generated by the functional ()5.6|) . Thus we 
find that the Hamiltonian functional is a conservation law of Eqs. (|5.7j) . Its 
continuous limit is the functional 7i. Physically, we interpret this conservation 
law as conservation of energy. 

2. The space translation X2 corresponds to the evolutionary operator (|4.7|) . 

which we already examined. This symmetry leads us to the conservation law 

V2 (see (gSl)- 
h 

3. The evolutionary symmetry X^ is generated by the functional 



2 

n 

which is a conservation of mass for the semidiscrete system (|5.7() . In the con- 
tinuous limit it goes into the functional 

V3 = ^ f {v^ + w^) dx. (5.9) 
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4. Galilean transformation X4 is not admitted in the semidiscrete case. Its 
coefficients violate the second mesh invariance condition in ()3.4() . The action of 
the generated by transformation on the independent variables 

t = t; X = X + 2ta 

clearly shows that it destroys the mesh geometry. It follows that the continuous 
conservation law corresponding to the movement of the center of mass 

P4 = y (^{v'^ + w"^) + tiwv^ - vw^)^ dx (5.10) 

has no counterpart in the semidiscrete case. 

5. The additional symmetry X5 is admitted by equations (|5.7|) with quadratic 
F, but it is not Hamiltonian (both in the continuous and discrete cases). 

We sup up the considered symmetries and their generating functionals in 
the Table 2. 



6 Conclusions 

In the paper we considered semidiscrete canonical Hamiltonian equations and 
showed how to find conservation laws of such equations using Noether's theo- 
rem. Our interest to canonical Hamiltonian equations is also motivated by their 
connection with Euler-Lagrange equations :3j;. 

Many PDEs can be presented as Euler-Lagrange equations for appropriate 
Lagrangian functionals J7] and can be rewritten as canonical Hamiltonian 
equations (ll.lj) . Let us illustrate this on the example of the nonlinear wave 
equation (|4.1() which is the Euler-Lagrange equation 

— - — - — D — D — (6 1) 

6v ' dv dv dvt ^ dvx 

for the Lagrangian functional 

C = 1 1 L{v,vt,v,)dxdt, L = ^-^-Viv). (6.2) 

With the help of the Legendge transformation 

dL f [ dL \ 

w = — = vt, H= [ — - L \ dx 
OVt J \dvt ) 

we obtain the Hamiltonian functional ()4.3|) . It generates equations ()4.2|1 , which 
are equivalent to (|4.1|1 . 

Similar connection between Euler-Lagrange equations and canonical Hamil- 
tonian equations can be established in the semidiscrete case. Let us consider 
the semidiscretization of functional 1)6. 2() 

^= / E^^^'^'^i)^'^*' L = ^-^-Viv). 

h J ^ h h h 2 2 
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Its Euler-Lagrange equation 



has the form 



5 d ^ d d 

dv ov ov ov -hOV I 

h 



v = V2-V'{v). (6.4) 

h 



We can introduce the discrete Hamiltonian functional 

dL 



w = -^ = v, = yZ{-^^ - l] h, 

ov ^ n ^ 

which in our case is the functional (|4.4|) . It generates the semidiscrete evolution 
equations ()4.5() . which are equivalent to H6.4|) . 

Thus, looking for conservation laws of the semidiscrete Euler-Lagrange 
equations, one can consider the equivalent canonical Hamiltonian systems and 
find conservation laws in the Hamiltonian framework. They have the form 

Vdx and Vp/i (6.5) 

in the continuous and semidiscrete cases correspondingly. If necessary, the con- 
servation laws can be rewritten in term of the original variables. For example, 
the conservation laws of the semidiscrete wave system (|4.5|) with an arbitrary 
potential V{y) found in Section 4 can be rewritten in terms of variables used 
in the Lagrangian approach 



dL 



For a quadratic potential we also have the infinite series of the conservation 
laws 

Rk = yZvDo{v2k)h = -y^vDo{v2k)h, k = 0,1,2,... 
n ^ ^ a ^ ^ 

and the conservation laws 

r = 7 {a(t,x)v — at(t,x)v)h, 

where the function a{t,x) must satisfy the equation (|4.im) . 

One of the advantages of Lagrangian approach over Hamiltonian one is a 
possibility to find local conservation. In the Hamiltonian framework we can 
find only global conservation laws of the form ()6.5() . This advantage of the La- 
grangian approach is difficult to represent for the semidiscrete equations because 
as we have seen some conservation laws have densities which involve discrete 
presentations of the continuous derivatives and, consequently, are not local. 
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For simplicity we considered only the case of one-dimensional space. The 
extension to the multi-dimensional space is straightforward. 
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Appendix 



Table 1. A number of conservation laws for the wave system (4.5) without a 
potential V'{v) = 0, which correspond to particular cases of the symmetry 



Function a(t, x) 


Operator Z 


Discrete conservation law 


Continuous conservation law 


1 


d 

dv 


n 


J wdx 


t 


t ^ ^ 
dv dw 


^^(iw — v)h 
n 


J {tw — v)dx 


X 


d 

dv 


xwh 

n 


Jxwdx 


tx 


d d 
tx—+x— 
dv dw 


'y^^{txw — xv)h 
n 


Jitxw-xv)dx 


t'^ + x^ 


/ 9 9^ d ^ d 


J2ii't'^ + - 2tv)h 
a 


J {{t'^ + x'^)w - 2tv)dx 
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Table 2. Some symmetries and their generating fiinctionals for the canonical 
bracket. The functionals (up to some time-dcpcndcnt functions) are 
conservation laws of the semidiscrete canonical Hamiltonian equations. 



Operator 


Operator in evolutionary form 


Functional 


Continuous limit of functional 


d 
di 


h \ dw j dv h \ dv j dw 


H 

h 


n 


d 

dx 


^ , ^ d ~ , . d 
Do{v)^+Do{w)^ 

h OV h OW 


wD o{v)h 


J wvxdx 




d d 
dv dw 




^ / 




h h dv h h dw 


y2wDQ{v2k)h 
n ^ ^ 


J wv^2k+i)dx 




a{t,x)-^ + at{t,x)-^ 


^^(a(t, x)w — at(f, x)v)h 
a 


Jiait,x)w-a,it,x)v)dx 
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